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ABSTRACT 

Whether a short gamma-ray burst (GRB) is caused by a black hole (BH) or a neutron star (NS) 
after the merger of a NS binary is a crucial problem. We propose a BH model that explains short 
GRBs with long-lasting activities such as extended emission and plateau emission up to 10000 s. To 
extract the BH rotational energy, the topological evolution of the magnetic field should accompany the 
mass ejection, mass fallback, and magnetic field reconnection. The observations suggest the magnetic 
held decay from ^ 10^'^ G to ^ 10^^ — 10^^ G at the BH, bounded below by the pre-merger strength 
and kept constant while the luminosity is constant, and the fallback mass of ^ 10“"'^ — 10 “^Mq, 
comparable to the ejecta mass implied by the macronova (or kilonova) in GRB 130603B. The BH 
model has implications for gravitational waves and the equation of state of NS matter. 

Subject headings: - 


1. INTRODUCTION 

A black hole (BH) or a neutron star (NS)? This is the 
fundamental question about the ce ntral engine of short 
gamma-ray bursts (GRBs; e.g., I Zhang fc Meszaros I 
12004 [Berger Il20lll . A leading model for the short GRBs 
is a binary NS merger. A merger produces a massive NS, 
which may or may not collapse to a BH depending on the 
equat ion of state of NS matter (|Bartos, Brady &: Markal 
I2013f) . In a BH-NS binary merger, the central en¬ 
gine would be the BH. Obviously, there is a signif¬ 
icant difference in the jet formation between a BH 
and NS. The central eng ine also affects the gra vita- 
tional wave (GW) signals (jHotokezaka et al.'1l2011ll and 
the e l ectromagnetic counterparts (e.g., Nalar^&PiranJ 
2014 RTao et al. Il2013l [Takami. Kvutoku fc loka 112014 

mi I [2 ' 


Kisaka. loka &: Takami IT2OI5II . Detection of these sig¬ 

nals will be soon realized by laser interferometers such 
as Advanced LIGO, Advanced VIRGO and KAGRA. 

X-ray observations suggest long-lasting activities of the 
central engines of short GRBs. Following the prompt 
GRB for ^ 10 ~^ — 1 s, some events sh ow extended emis¬ 
sion (~ 10^ s dBarthelmv et ;il~ll2005f) . which lasts much 
longer than the typical accretionEj. Furthermore, some of 
events accompany even a longer activi ty ("^ 10^ — 10^ s; 
hereafter we call this plateau emission) (Gomnertz et al. I 
I2013I iGomoertz. O’Brien fc Wvnn 11201^ FI The sharp 
drop of the light curve is produced only by the ac- 
tiyity of the central engine (jloka. Kobayashi fc Zhangl 
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^ We define the extended emission as the emission for ~ 
10^ s, which also inclu des the plateau com ponent analyzed by 
IRowlinson et al~l l|2013l j and ILii et ak~l 1120151) . 

^ The plateau c omponents with long; tim escale (~ 10^ — 10^ s) are 
also identified bv IRowlinson et al. I 1)20131 ). Note that the plateau 
emission would be sometimes hidden by the afterglow emission or 
below the detection limit. These eve nts may correspond to the “no 
breaks” in IRo wlinson et al. I 1)20131) and “no plateau samples” in 
ILii et aL~l 1)20151 ). 


2005lj. In this Letter, we consid er the samples in 


Gompertz. O’Brien fc Wynn I (12014) . 


To explain the long-lasting actiyities, a highly magne¬ 
tized and rapidly spin ning NS model (magnetar model) 
is co nsidered (e.g., iMetzsrer. Quataert fc Thompson] 
HOOI). This comes from the fact that the spin-down 
timescale is ~ 10^ s for a NS with the dipole magnetic 
held ^ 10^® G and the initial rotation period ~ 1 ms. 
The spin-down luminosity is also comparable to the ob- 
seryed one within a reasonable range of the parameters 
0 . 


In the case of the BH engine, fallback accretion is 
considered to produce the long-lasting activities (e.g., 
iLee fc Ramirez-Ruiz Il2007^ . In fact, the released energy 
is enough to power the observed luminosity. However, 
the accretion rate of the fall back matter foll ows a single 
power-law (see Equation (g]); lRosswogl2007D . Then, the 
resulting activity does not seem to explain the character¬ 
istic timescales in the light curve. 

In this Letter, we present a BH model that repro¬ 
duces the long-lasting activities of short GRBs. The 
point is the topological evolution of the magnetic 
held associated with the mass ejection and fallback 
processes. We show that the mass ejection, mass 
fallback, and magnetic held decay via reconnection are 
all inevitable events fo r the jet formation by Blandford- 
Znajek (BZ) process ([Blandford fc Znaiek 1 119771 ). We 
hnd that the magnetic held changes to ~ 10^"^ G, 
^ 10^^ G, and ^ 10^^ G in the prompt, extended and 
plateau emission, respectively. The required mass of 
the fallback matter is ^ — 10 “^Mq, which is 

consistent with th e recent n umerical simulations (e.g., 
Hotokezaka et alH 120131: iKvutoku. loka fc Shibatcil 
20131 : iKvutoku et al. I 1201^ as well as the observed 
infrared excess (the so-called macronova 0 or kilo- 


® The prompt GRB jet should be launched before the outflow of 
the long-lasting activities because the outflow becomes too thick 
for the jet to propagat e keeping the prompt timescale 0.1 s 
( NagakuTa_et_an |20ij) unless the outflow is u nreasonably cold 
I EzzoUa &; Kumar II2015I : ICiolfi &; Siegel II2015I) . 

® We use the term “macronova” as a transient with a binary 
NS merger, especially thermal radiation from the merger ejecta. 
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nova) associated w i th short GRB 130603B (e.g. , 
Tanvir et al. I [2M^ i B erger. Fo ng fc Ch ornock I 120131: 
T akami. Nozawa &: loka 112014 iKisaka. loka &: Takamil 
201^ . We introduce our BH model in Section 2. In 
Section 3, we compare the theoretical light curves with 
observations. Discussions follow in Section 4. 


2. BLACK HOLE MODEL 

In our BH model, a BH is the central engine of 
short GRBs including long-lasting activities such as ex¬ 
tended emission and plateau emission. We consider 
the BZ process for the BH to launch a relativistic jet 
(jBlandford fc Znaiek 1119771) since the neutrino-neutrino 
annihilation is not effective at late time (> 1 — 10 s). 
Our model is based on general topological consideration 
in Figure [U without resorting to specific processes such 
as radiation mechanisms. 

There are three key ingredients for the BZ process, (i) 
rotation of the BH, (ii) a magnetic field strength on the 
BH, and (iii) large-scale, poloidal configuration of the 
magnetic field, which means that the characteristic size 
of the poloidal field on the BH is much greater than th e 
outer light cylinder (jBeckwith. Hawlev fc KroliFll2008D . 
We consider the BH with mass Mbh, a spin parame¬ 
ter a = Jc/GMbh and a magnetic flux 4 >bh 
where J is the angular momentum of the BH, c is 
the light speed, G is the gravitational constant, Hh 
is the strength of the magnetic field at the BH and 
rn is the radius of the BH horizon. Then, the total 
power of the BZ jet is (e.g., [Blandford fc Znaiek Ill977t 

[Tchekhovskov. Naravan fc McKhmev 1120111) 


Lbz ' 


47rc 




H^BHi 


where k 0.05, the angular frequency of the BH is 

a* *c 


Dh = 


2rH’ 


( 1 ) 


( 2 ) 


and a, = a/M bh is the dimensionless spin parameter. 

Phases I - HI in Figured] The rotational energy of the 
BH formed after the binary NS merger is huge. After the 
inspiral phase of the binary NS merger (phase I), a hyper- 
massive neutron star (HMNS) is formed (phase H), whose 
gravitational collapse is prevent ed by differential rotation 
and thermal pressure (e.g., iBartos, Bradv fc Marka I 
I2013D . Within the transport timescales of angular 
momentum and thermal energy (< 10“^ — Is), the 
HMNS eventually collaps es to a BH with its surrounding 
torus (phase HI) (e.g ., iBartos, Bradv fc Marka 1 120131 
iSekiguchi et al. 1120151) . From the numerical simulations, 
the dime nsionless spin parameter of t he collapsed BH is 
a* ~ 0.7 (|Shibata fc Taniguchi 11200^ . Then, the avail¬ 
able rotational energy of the BH is 


Ernt. — I 1 — 


1 + — al 


Mbbc^ 2 X lO'^^ 


which is enough to explain the total energy of a short 
GRB. The problem is how to extract the rotational en¬ 
ergy up to long timescale, ~ 10“^ s. 


inclu ding an engine-powered macronova IKisaka, loka fc 'I'akamfl 

l20iB . 


A magnetic field is amplified after the merger. In pre¬ 
merger phase, a non-recycled NS has a dipole magnetic 
field i?NS ~ 10^^ G (phase I). Since the radius of NS, 
i?NS ~ 10 km, is almost equal to that of the BH horizon 
with the mass Mbh 3Mq, the strength of the magnetic 
field at the BH is about Bb ^ Hns thanks to the flux 
conservation. The corresponding BZ power is too weak 
to explain the typical observed luminosity of the prompt 
emission of short GRBs, ^ 10“^® — 10®^ erg s“^. At least 
the strength of the magnetic field Bb ^ 10^'^ G is re¬ 
quired for the BZ process in Equation O- Therefore, the 
magnetic field should be amplified by some mechanism 
such as Kelvin-Helmholtz instabilit y and/or magnetoro- 
tational instability (phase H; e.g., IShibata et al. l[2Mll : 
iKiuchi et al. 112014) in matter such as the HMNS and/or 
the torus. Re cent numerical simu lations support such a 
picture (e.g., IKiuchi et al. 11201^ . 

The magnetic field should be expanded to large 
scale for the BZ process to work. At the same time, 
this expansion should be associated with the mass 
ejection because of the frozen-in condition (phase HI). 
The mass ejection is caused by suc h as winds driven 
by d yn amical interactions (e.g., iHotokezaka et aT~1 
120131 : iKvutok u . loka fc Shibata I 120141 1. neutrinos 
(e.g., iDessyt et al. I 20091). magnetic field (e.g. . 


IShibata et al. II2011I : iPaschalidi s. Ruiz fc Shapi ro 1201'^ 

and yiscous heating (e.g., iFernandez et al. I 1201511 . 
Although these scenarios do not currently succeed 
in numerical simulations, the amplification and the 
expansion of the magnetic field are necessary to explain 
the obseryed short GRBs regardless of details. 

Phases IV - VH in FigurfJT] The central BH accretes 
mass from the surrounding torus. The timescale of the 
accretion is determined by the yiscosity timescale at the 
inner radius of the torus Gis 0.1 s (phase IV; e.g., 
iZhang fc Meszarosi 1200411 . It has been widely discussed 
that the torus accretion could explain the prompt emis¬ 
sion of short GRBs. 

In addition, the ejected matter falls back to the BH. 
This is because the ejecta beco me homologous due to the 
interaction within ejecta (e.g., IHotokezaka et al.~ll2013ll 
and the interior part of the ejecta does not exceed the es¬ 
cape velocity (phase V). The fallback mass is comparable 
to or larger than that of the escaping ejecta because the 
inner part of ejecta is more massive (jHotokezaka et al. I 
120131) . The temporal evolution of the mass accretion rate 
is described by the following power-law function (e.g., 
IRosswo^|200'A . 


M = 


2Mf / t 

3 tvis V tvi 


-5/3 


{t > tvis)5 


(4) 


where Mf = Mdt is the total fallback mass. 

*vis 

The fallback enables the long-lasting activities of the 
BZ jet because the pressure of the fallback matter sup¬ 
ports the magnetic flux on the BH, (phases V and VH). 
As long as the magnetic flux 'Fbh is constant, the BZ 
power does not directly depend on the mass accretion 
rate (in Equat ion]!]) and the light curve becomes plateau¬ 
like Lbz ~ (|Tchekhovskov fc Giannios II2015D . 

The important point is that the fallback matter drags 
the magnetic field line to the BH because of the frozen- 
in condition and eventually forces the anti-directed mag- 
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Fig. 1.— Schematic pictures of our BH model for short GRBs. See Section 2 for details. 


netic fields to reconnect (phases IV and VI). Then, the 
magnetic flux on the BH Tbh as well as the BZ power 
Lbz decrease. In other words, the long-lasting activities 
require the fallback matter, which itself inevitably leads 
to the magnetic flux decay. This topological argument 
does not depend on the specific physical processes. Dif¬ 
ferent values of the magnetic fields explain the BZ power 
at different phases, prompt (phase III), extended (phase 
V) and plateau emission (phase VII) as quantitatively 
shown in the next section. The minimum magnetic flux 


after the reconnection (phase VII) is determined by the 
initial flux of a NS before the merger (phase I). 

Phase VIII in Figure [TJ The BH activity ends if the 
pressure of the fallback matter becomes too small to sup¬ 
port the magnetic flux. Then, the magnetic field lines 
escape from the BH (phase VHI). 

3. LIGHT CURVE 

Figure [2] shows the theoretical light curve with a thick 
curve. Black lines denote the constant luminosity phases 
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Fig. 2.— A representative light curve for prompt, extended 
and plateau emission in our BH model. Observational data of 
GRB 070714B is obtained from UK Swift Science Data Cen¬ 
tre. Time shown in the horizontal axis denotes the rest-frame 
time since Swift/BAT t riggers. For the redshift value, we follow 
ICompertz et alTl II2Q13I) . The number III - VIII corresponds to 
the phase in Figure^ For the rebrightening component at 1 s, 
we consider flaring activities discussed in Section 4. 

(Ill, V and VII), and light-green lines denote the decay 
phases (IV, VI and VIII). Constant luminosities at three 
phases (ill, V and VII) suggest constant magnetic fields 
in Equation O; Bn ^ 10^"^ G at the prompt emission 
(phase III), ^ 10^^ G at the extended emission (phase V) 
and ~ lO^^G at the plateau emissi on (phase VII), respec- 
tively. Here, we assume a* ~ 0.7 (IShibata &: Taniguchi I 

l200fili . so that TH = (1/2) ^I -|- Rs ^ 0.86i?s, 

where i?s = 2 GMbh/c^ is the Schwarzschild radius. To 
convert the luminosity from the BZ power Lbz to the 
observed isotropic luminosity L, we take into account 
the beaming correction (d? ~ 10~^: l^ng et al.~ll2014fl 
and t he radiative efficiency (ry ^ 0.1; e.g., iZhang et al. I 

l200l . 


L^r]{2/ef)Lnz-W^Lnz- (5) 


Using Equations o and ([5]) , the strength of the magnetic 
field is determined by the observed luminosity L as 


Hh~3 X 10^ 


'v/df 

102 


- 1/2 


/ Mbh \ 


-1 


L 


1/2 


\ lO^^erg 


,-i 


G. 


( 6 ) 


A characteristic timescale of the BZ jet activ¬ 
ity is determined by the pressure balance near the 
BH. As the mass accretion rate decreases, the pres¬ 
sure of the fallback matter pf falls short of the 
magnetic press ure pb (making a magnetically-a r rested 
disk; e.g., | Bisn ovatvi-Kogan fc Ruzma ikin I [197^ 
iNaravan. Igumenshchev &: Abramowicz I I^OSfb Then, 
the magnetic flux expands the torus and decreases on 
the BH, leading to the reduction of the BZ power. The 
pressure of the magnetic field is 


f —1 

8tt \rnj 


(7) 


under the magnetic flux conservation, where R is the 
radial distance. The pressure of the fallback matter is 


GMnn M 
2TrR^VB, 


( 8 ) 


where vji is the radial velocity. For the radial veloc¬ 
ity 'f^R) we assume ur = evs where vs = 
is the free-fall velocity. For the value of e, we 
adopt e ^ 10“2 which is supported by the obser¬ 
vations an d numerical simulations of the relativistic 
jets (e.g., [T chekhovs kov. Naravan fc McKinney 1 120111 : 
iZamaninasab et al. [1201411 . Then, the magnetospheric 
radius Rm, where the equilibrium point between two 
pressures pi (Equation [7]) and pn (Equation [8]) , is 



Using Equation m and the temporal evolution of the 
mass accretion rate (Equation 0]) , we obtain the charac¬ 
teristic timescale of the BZ jet 



Using Equations and m, the total fallback mass is 
determined by the observed luminosity L and duration 
T as 


Mf ~ 1 X 10"2M( 


® ( 10 - 2 ) 


'??/6»2 

102 


X 


\ 10"^^erg s“i 




-2/3 


( 11 ) 


After the time t > T, the BZ power evolves as Lbz oc 
/-40/9 derived by the magnetic flux ^^bh oc i?“2 and the 
time dependence oc from Equations (0]) and dH]). 
Thus, we model the BZ power as 


Lbz (X + . (12) 


The maximum value of the BZ power Lbz is deter¬ 
mined by the accretion power Mc^. In fact, using the 
condition pf > pn and Equations o, © and ©, the 
ratio is Lnz/{Mc^) < 1. For comparison, we plot 
the mass accretion rate with the beaming correction 
(2/02 )Mc2 10^Mc 2 in Figures 0] and [3] as thin dot¬ 

ted lines. 

Our theoretical light curve is consistent with the obser¬ 
vations. As an example, we plot the observational datcQ 
of short GRB 070714B, which has extended and plateau 
emission in Figure [3J 


http://www.swift.ac.uk/index.php 
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Fig. 3.— Theo retical light curves for 9 s hort GRBs used in lGompertz et al. I ll2Q1 3l):IGomDertz. Q’Brien &: Wynin 11201411 . For the redshift 
values, we follow [Gompertz et al. I ||2013I'1 except for GRB 051227 iz = 0.8: ID’Avanzo et al"Tl200^ 


In Figure [3l our model is consistent with other 9 short 
GRBs used in lGompertz. O’Brien fc Wynn I (|2014f) . The 
ratio Lsz/(Mc^) ^ 1 is satisfied for all samples. 

4. DISCUSSION 

Figure 0] shows the total fallback mass Mf and mag¬ 
netic helds Rh at extended (blue triangles) and plateau 
emission (red circles) evaluated by the observed luminos¬ 
ity L and duration T with Equations (O and (fTTl) . The 
magnetic field at the plateau emission is Rh ^ 10 ^^ —10^^ 
G, which is consistent with a typical value of pulsars 
and the observed value of a non-recycled pulsar, PSR 
J073 7-I-3039B, in the double pulsar system (|Lvne et al. I 
l2004fl . For the extended emission, slightly strong mag¬ 
netic field Rh 10^^ G is required to explain the ob¬ 
servations. Different values of the magnetic field be¬ 
tween prompt and extended emission may suggest that 
the magnetic field in the torus formed at the collapse of 
the HMNS is larger than that in the ejecta. On the other 
hand, the total fallback masses Mf for both extended 


and plateau emission in Figure |4] are consistent with the 
ejecta mass Mf ^ 10 ~"^ — 10~^Mr7^ obtained from the 
numerical simulations (iHotokezaka et ai~1l2013[l and the 
observation of a mac ro nova following GRB 130603 B (e.g. , 
Tanvir et al. I I2013I i B erger. Fong &: Ch ornock I [20131 
Takami. Nozawa fc loka I l2014t iKisaka. loka fc: Takami I 
201511 . The fallback masses for the extended and plateau 
emission are similar in each event, partly supporting our 
model. 

In Figure 4, we also compare our BH model with GRB 
080503, which has a bright extended emission and a 
weak or no plateau emission. Th e parameters were unob¬ 
tainable in the magnetar model (iGompertz et al. 1120131 : 
iGompertz. O’Brien fc Wynn 11201411 . We obtain only the 
upper limit on the magnetic field Rr at the plateau phase 
(red arrow), which are consistent with the other events. 

Rapid declines of some light curves may indicate that 
the magnetic held is decayed by reconnection since Equa¬ 
tion (021) assumes the flux conservation. Note that the 
released energy due to the reconnection is negligible for 
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Fig. 4. — Estimated total mass of the fallback matter Mf 
and magnetic field Bn for extended emission (blue triangles) and 
plateau emission (red circles). Solid lines connect each event of 
short GRBs. 


the extracted energy by the BZ process. 

Some events are ac companied by X-ray flares (e.g., 
iMargutti et al. I l2011h . These activities could be ex¬ 
plained by the accretion of blob with the same mag¬ 
netic pol arity as d/RH. the decr ease of the jet opening- 
angle 9j (iMizuta fc loka I [201311 and/or the increase of 
the radiative efficiency rj. After the time t > T, the 
accretion could be episodic, so tha t the magnetic flux 
on th e BH 'f'sH fluctuates (e.g., iProga fc Begelmanl 
2003f). This potential ly gives the flaring activities (e.g., 
Proga fc Zhang 1 12006(1 . 


The extended or plateau emission is only seen in a 
fraction of the short GRBs. As mentioned in Section 1, 
the plateau emission could be hidden by the afterglow 
emission or below the detection limit. Outflows from 
the acc retion disk could also i nteract with the fallback 
matter (jPernandez et al. II2015I1 and reduce the extended 
and plateau emissions. 

The model in Figure [T] is applicable to the merger of 
a NS-BH binary. The merger produces a BH-torus sys¬ 
tem with mass ejection (e.g ., [K yutoku . loka &: S hibat'^ 
I2013I iFoucart et al. 112015 : iKvutoku et al. Illoislh Dif- 
ference from the case of the binary NS merger is that a 
HMNS is not formed in the NS-BH merger case. How¬ 
ever, the m agnetic field amplification in the tor us is pos¬ 
sible (e.g., iPaschalidis. Ruiz fc Shapirol 1201411 . There¬ 
fore, the phases HI - VHI would be the same for both 
cases. 

The model in Figured] is also applicable to long GRBs 
since the mass ejection and fallback could also occur with 
the supernova explosion and/or the central engine activ¬ 
ities. This topic will be discussed in future work. 

The long-l asting activities would significa. ntlv affect the 
macronovae (|Kisaka. loka fc Takami 11201511 since the en¬ 
ergy injected to the preceding ejecta suffers from rela¬ 
tively small adiabatic cooling. This issue will be studied 
in a separated paper. 
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